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1.0 INTRODUCTION 


Encounter with wind shear during approach and takeoff has gained 
recognition within the past few years as a major hazard to flight safety. 
Several accidents and near accidents have now been officially attributed 
to wind shear. There were probably many more prior to instruments and 
basic knowledge of the atmosphere which made identification of wind 
shear as a contributing factor possible. These accidents and near 
accidents have been mistakenly resolved as pilot error. The most notable 
accident with wind shear as a major cause involved the Eastern 66, 

Boeing 727, which crashed on approach in a thunderstorm wind environment 
at John F. Kennedy International Airport on June 24, 1975, resulting in 
113 fatalities [1].* A more recent incident took place on August 22, 
1979, in Atlanta as Eastern 693, also a Boeing 727, was forced to abort 
approach and initiate a go-around after a relatively severe loss of 
altitude and airspeed [2]. 

Wind shear is a spatial or temporal gradient in wind speed and/or 
direction. It has many sources in the atmosphere. A common but 
generally benign example is the logarithmic wind profile associated with 
boundary layer flow over homogeneous terrain. A vertical shear exists 
in the longitudinal wind--the wind velocity has some magnitude w at 

X 

height h which decreases logarithmically to zero at the ground. Other 
sources of shears are cold and warm fronts, thunderstorms, and perturbed 
flows around obstacles. 

To illustrate the effects of wind shear on an aircraft, two simple 
examples are examined. Since an airplane's flight depends on horizontal 
airflow (relative airspeed) over its wings to generate lift, any changes 
in the relative airspeed will alter the balance of forces that keep the 
aircraft in straight, unaccelerated, level flight. Therefore, a tail 
wind shearing to a head wind will increase the relative airflow over the 


*Numbers in brackets correspond to similarly numbered references in 
the List of References. 


wings, increasing the lift force, and causing the aircraft to ride high 
on the glide slope. In the case of a head wind shearing to a tail wind, 
the relative airflow over the wings will decrease, the lift will decrease, 
and the aircraft will fall below the glide slope. This is a particularly 
hazardous condition since reaction time and engine spool -up time become 
relatively small as the plane approaches the runway threshold; premature 
ground impact may result. Escape from this hazardous situation depends 
on the severity of the wind conditions, the physical and performance 
characteristics of the aircraft, and the skill, experience, and reaction 
time of the pilot. If a trained pilot recognizes his decrease in air- 
speed and deviation below the glide slope, he can apply the necessary 
thrust and elevator controls to perform a go-around or land the plane 
safely. 
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2.0 WIND SHEAR STUDIES 


Previous computer simulations of aircraft flight through wind shear 
have suggested that the phugoidal oscillations of commercial jetliners 
are excited when the wind spectral frequency approaches the phugoidal 
frequency. 

McCarthy, Blick, and Bensch [3] used aircraft measured thunderstorm 
environment winds as input to their numerical aircraft flight simulation 
model. They used a two-dimensional linear system of aircraft motion 
equations to perform fixed-stick (constant thrust setting and elevator 
deflection angle) flight simulations with an aircraft characteristic of 
a Boeing 727. The results of the simulations show that horizontal gusts 
at approximately the phugoidal frequency are responsible for most of 
the flight path deviation and airspeed oscillation. It was also noted 
that flight simulations through a wave-form longitudinal wind profile 
representative of the winds encountered by Eastern 66 [1] resulted in 
poor approaches. 

Frost and Crosby [4] numerically simulated fixed-stick and auto- 
matic controlled flight of aircraft characteristic of a DC 8, Boeing 
747, DHC-6 Otter, and STOL (short takeoff and landing) through thunder- 
storm wind fields. They used a nonlinear two-dimensional three-degrees- 
of-freedom model. Two-dimensional thunderstorm wind field data were 
provided in tabulated format with a computer lookup routine. The 
longitudinal and vertical wind component data used in the study were 
measured with an instrumented 500 m tower at the National Severe Storms 
Laboratory at Norman, Oklahoma. Time history measurements were converted 
to a horizontal spatial coordinate grid using Taylor's hypothesis for 
use in the computer lookup routine. Frost and Crosby's simulations 
revealed highly amplified phugoidal oscillations along flight trajectories 
through these storms. 

Extensive work was done in simulated automatic controlled flight 
through boundary layer shear, perturbed wind fields around buildings, 
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and thunderstorm gust fronts by Frost and Reddy [5]. A feedback control 
loop was developed which passes airspeed error, horizontal and vertical 
acceleration, and pitch command signals through variable gains to drive 
the throttle and elevator controls. A detailed discussion of this auto- 
matic control system appears in Appendix A. 

The aforementioned work has established a relationship between air- 
craft response to longitudinal thunderstorm gusts; however, the fact 
remains that influence of a human pilot has not been included. The 
fixed-stick method is not likely to be flown in a thunderstorm environ- 
ment, whereas the automatic control system is too perfect. To achieve a 
more realistic simulation, this study examines the response of the 
pilot-aircraft system to hazardous wind environments. 

The work of Frost and Reddy [5] is expanded to include a human 
pilot as controller. Computer-simulated flight through variable winds 
measured along typical flight paths is carried out. The wind data were 
measured along the glide slope with a Doppler radar. Simultaneously, a 
Queen Air aircraft made controlled approaches along the glide slope. 

The measurements were carried out as part of the SESAME '79 program, 

April and June, 1979, in Oklahoma. 

The overall objective of the flight experiment program is to test 
the feasibility of a Doppler radar wind shear warning and detection 
system. The operational system concept is to aim a 10.2 cm pulsed 
Doppler radar along the glide slope at commercial airports and thus 
supply air traffic control personnel with real-time wind information on 
the approach path. The Doppler radar system is capable of making a 
complete scan along the glide slope from 3.5 km (minimum range) to 22.5 
km (maximum range) and obtain almost instantaneous wind information in 
150 m steps. Such a system would be extremely beneficial in detecting 
hazardous shear conditions for any aircraft landing on a monitored 
approach path. Air traffic control personnel can monitor worsening wind 
shear conditions between takeoffs and landings of various aircraft. 

Also, a real-time spectral analysis can be performed on the Doppler wind 
data to determine if phugoidal frequency shear waves particularly hazard- 
ous to a given aircraft type exist. The phugoid oscillation is discussed 
in Section 3.0. 
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In order to compare the effects of the Doppler-detected wind fields 
on aircraft other than a small commuter-type Queen Air aircraft, an 
aircraft having characteristics of a Boeing 727 was also flown through 
one of the wind fields. In addition, the aircraft was flown through 
sinusoidal longitudinal wind waves having a frequency near the phugoid 
frequency of the aircraft. The wind field was constructed to have a 
head wind shearing to a tail wind, thus giving an effect similar to 
flying directly through the downburst region of a thunderstorm cell. 
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3.0 THE PHllGOID OSCILLATION 


An aircraft responds to a longitudinal disturbance with oscilla- 
tions in speed, vertical height, and pitch rate. The oscillatory motion 
can be broken down into two components: the short- period oscillation 

and the long-period phugoid oscillation. The short-period oscillation 
is heavily damped, and its associated pitch rate changes "wash out" 
shortly after the disturbance is encountered. The aircraft is then left 
with a long-period, low-damped phugoid oscillation which produces slow 
changes of speed and height. This is usually not a serious problem for 
a pilot to control since the speed and height changes are not severe and 
the time period is long. It is insidious, however, when the pilot's 
attention is directed to other factors. Moreover, encounter with a 
large wavelength longitudinal disturbance (horizontal wind shear) may 
exaggerate the phugoid mode, causing a control problem for the pilot. 

Analysis of the phugoid oscillation assumes small pitch' rates, 
constant angle of attack, negligible drag effects, and a constant sum of 
potential and kinetic energy. A simple example involves an aircraft 
encountering a head wind gust. The initial response is an increase in 
lift, a pitch-up of the nose, and an increase in height; kinetic energy 
is converted to potential energy, causing the airspeed to drop and the 
nose to pitch down. Potential energy is converted back to kinetic 
energy as the plane begins to descend, picking up airspeed until it once 
more begins to climb. The kinetic energy of descent is converted to 
potential energy of height. This simple harmonic motion will damp out 
after some time. 

From Dickinson [6] the degree of phugoid damping for a particular 
aircraft depends on the aerodynamic derivative x , which is the rate of 
change of force in the x-axis (direction of flight) due to changes in 
forward speed. This is related to the drag characteristics of the 
aircraft. Also known as the speed damping derivative, x u is usually 
negative and proportional to -Cp where is the drag 
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coefficient. If the aircraft is disturbed in such a way as to increase 
its forward speed, lift will increase as well as drag, which will eventu- 
ally reduce the speed back to its initial value. 

From McCarthy, Blick, and Bensch [3] the exponential phugoid damping 
period, or the time for the oscillation to decay to 1/e of its initial 
amplitude, is given by T = Va 0 /g(C|_ 0 /C[) 0 ) where V aQ is equilibrium 
airspeed (m/s), g is the gravitational acceleration (m/s 2 ), and C|_ Q and 
Cq q are the lift and drag coefficients for the particular 
aircraft. Thus, the high approach speed (V aQ ) and high lift-to-drag 
ratios of commercial jet transport aircraft contribute to a large expo- 
nential phugoid damping period. This simply means that the phugoid 
oscillation of a jet transport will not damp out as quickly as that of a 
smaller aircraft. This adds significance to the observation that smaller 
aircraft have escaped longitudinal gusts that later caused large commer- 
cial aircraft to crash. 
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4.0 THE PILOT MODEL 


The pilot models used in the computer-simulated flights are 
described in this section. 

For the preliminary study, a parametric study of pilot performance 
rating was carried out. The pilot performance rating model was incor- 
porated into the control loop of Frost and Reddy [5] in the following 
manner. The automatic feedback control system as described in Appendix 
A evaluates the difference between the actual and desired state of the 
simulated aircraft on approach. The system monitors airspeed, horizontal 
ground speed, vertical velocity, flight path angle deviation and displace- 
ment from the glide slope and processes these parameters through variable 
gains to determine the necessary signals to the throttle and elevator 
servos. The servos adjust thrust and elevator angle to bring the aircraft 
back to the desired state. A pilot, who is assumed to be monitoring 
these same parameters at the same rate as the autopilot, would not be 
able to respond as fast as the autopilot to the information he receives. 
Thus, the pilot would not move the servos to correct the plane's devia- 
tions as efficiently as the autopilot. In effect, the pilot is always 
lagging the autopilot. Thus, the pilot's control signal inputs are 
reduced by a performance rating or "perfection percentage" gain where 


Initial 
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0 percent corresponds to zero inputs (fixed-stick) and 100 percent 
corresponds to "perfect" autopilot control. This is accomplished in the 
program by use of the difference equation 

y n = Vl - (y n-l - x n )k 

where y R and y p -j are the control signal outputs for the present and 
previous time steps, respectively, x p is the input control signal for 
the present time step, and k is the "perfection percentage" gain. The 
pilot constantly attempts to return the aircraft to the desired state, 
but this occurs at a slower response rate than the "perfect" automatic 
control system. 

•Classifying pilot response by means of a performance rating encom- 
passes the many intangibles encountered in pilot modeling which are too 
complex to simulate. These intangibles include pilot personality, 
training, knowledge and warning of the encountered wind shear, as well 
as the element of surprise. Hence, a pilot with a low performance 
rating (for example, 0.03, which corresponds to a minimal control input) 
may be classified as poorly trained, slow-to-react, unknowledgeable, or 
uninformed of the imminent wind shear. Clearly, more work needs to be 
done on pilot modeling specifically to determine pilot response to the 
wind shear environment. However for purposes of this study, the proposed 
pilot's "perfection percentage" gain gives useful results. 

In addition, two runs were carried out using human pilot transfer 
function data taken during simulator experiments by Adams and Bergeron 
[7]. They tested six pilots (ages 30 through 47) and two test engineers 
in a flight simulator equipped with an oscilloscope and control stick. 

Pilot workload was created by filtering computer-generated noise to 
obtain a simulated flight disturbance. While viewing the disturbances 
on the oscilloscope, the subjects were instructed to keep the oscillo- 
scope signal aligned with a marked reference point by moving the control 
stick. Adams and Bergeron measured the subjects' static gains, lead and 
lag time constants throughout the runs, examining variations between the 
subjects for given controlled dynamics (degree of vehicle controllability). 
The general form of the human pilot transfer function is 
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6 _ 

D 


k l T 


1 + — s 

T 


(t + s)2 


where 

6 = analog pilot output produced by moving the control stick (volts) 

D = input disturbance sensed visually by the pilot (volts) 
t = lag time constant (sec - ^) 

k-|/x = static gain (volts/volts), k^ units are volts/(volts)(sec) 
k 2 /x = lead time constant (sec), k 2 is dimensionless 
s = Laplace operator (sec - ^) 

This transfer function is useful for describing a pilot in a linear 
single-axis attitude control system in a compensatory tracking task 
(i.e., glide slope tracking) and provides a representative model for 
analysis purposes. 

Throughout the runs Adams and Bergeron [7] noted that k-| , k 2 , and t 
stabilized within 30 seconds to near constant values for the remainder 
of the run. Thus, k^ , k 2 , and x form a constant gains transfer function 
representative of a given pilot for a given run with specified controlled 
dynamics. A sample transfer function for pilot F (Table 1), assuming 
control dynamics = 1 (easily maneuverable and controllable aircraft), 
with k-| = 3 volts/(volts) (sec) , k 2 = 1 , and x = 4 sec - ^ is 


6 0. 75[1 + 0.25(s)l 

D LI + 0.25(s)J2 

The complete pilot-aircraft system response to a disturbance (in 
this case wind shear) is obtained by multiplying the pilot transfer 
function by the controlled system transfer function. For example, from 
Adams and Bergeron [7], multiplying the pilot transfer function 


6 
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and control system dynamics 
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TABLE 1. Transfer Function Data of the Eight Pilots Investigated in 
Reference [7]. 


Controlled 

System 

Dynamics 

Pilot 

Identifier 

k l 

volts 

(volts) (sec) 

6 T 

sec“T 

k 2 

(dimensionless) 

1 

A 

8.0 

8.0 

0.0 

1 

B 

6.5 

7.0 

0.0 

1 

C 

9.0 

11.0 

0.0 

1 

D 

5.0 

5.5 

0.5 

1 

E 

9.0 

10.0 

0.0 

1 

F 

3.0 

4.0 

1.0 

1 

G 

5.5 

6.0 

0.5 

1 

H 

3.0 

3.0 

1.0 


where 0 is the aircraft system output response (volts) to pilot output 
which is input to the system yields 


0 6 0 

D ~ D * 6 


k l T 

r k 2 

1 + — s 

t 


J + S) 2 


The general form of the pilot-aircraft transfer function 


0 _ 

D 


k r 

f k 7 

1 + — s 

T 

i 

+ S) 2 


can be rewritten as 


D 


k l k 2 




IT + s 
k 2 


+ s) : 


A transfer function of this form multiplied by a step response, 1/s, can 
be z-transformed into a difference equation of the form 

■ < C1 >Vl + < C 2 >V 2 + < D1 >Vl + < DZ >Vl 


(Neuman and Foster [8]) where y n is the output signal, y n ^ and y n ^ are 
the signal outputs of the last two previous time steps, respectively. 
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and x n ^ and x n ^ are the signal inputs of the last two previous time 
steps, respectively. The constants for this difference equation (Neuman 
and Foster [8]) are 


Cl = 2e 


■xt 


C2 = -e" 2rt 
D1 = — + e" Tt 

T 


(k,k 2 - 


k l !t - T 


and 


D2 = 


” (e _Tt - 1) - (k, k 2 - k ] )t 


-xt 


The response time t is set equal to 0.01 corresponding to the sampling 
interval of the computer program used in the flight simulations. 
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5.0 RESULTS AND DISCUSSION 


Doppler radar-detected glide slope longitudinal wind data for the 
twelve Queen Air SESAME approaches supplied by McCarthy and Elmore of 
the University of Oklahoma consisted of three sets of data for each 
flight. Two sets of data were Eulerian or "frozen" glide slope wind 
scans taken by the Doppler radar before and after the approach of the 
Queen Air down the glide slope. The third set of data was obtained 
with the ground-based Doppler radar sampling the wind just ahead of the 
aircraft as it flew down the glide slope. 

In the following simulations, only the Eulerian data sets are used. 
The raw Eulerian Doppler wind data were recorded in terms of slant range 
from the radar location and were then converted by McCarthy and Elmore 
to time history data by Taylor's hypothesis using the average approach 
speed of the Queen Air. Thus, 24 cases of time history wind data are 
available as input to the nonlinear flight simulation program. These 
cases were initially run in the fixed-stick mode to determine which 
cases, if any, had any significant wind shear to cause airspeed and 
flight path deviations. Generally, the wind shear in most cases was not 
severe enough to cause major deviations, but two cases revealed 
phugoidal-type oscillations in the aircraft's approach trajectory. One 
of these cases is selected for detailed examinations where the original 
fixed-stick trajectory is compared with autopilot- and pilot-controlled 
trajectories. In addition, the simulated flights of a B727 flown 
through sinusoidal phugoid frequency longitudinal wind waves of varying 
amplitude are presented. 

5 . 1 Simulated Queen Air Flight Through Queen Air Flight 17/Run 2 Wind 

Field 

The longitudinal winds measured during the SESAME Queen Air Flight 
17/Run 2 are used as input to the flight simulation program. Figure 1 
shows the head winds encountered during the flight as a function of 
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horizontal distance traveled. The initial increase and subsequent drop 
in head wind appears as a wave-form similar to those observed by 
McCarthy, Blick, and Bensch [3] and Frost and Crosby [4] which cause 
airspeed and altitude deviations during approach with fixed stick. 

Figure 2 shows the simulated fixed-stick trajectory of an aircraft 
having the characteristics of a Queen Air with the phugoidal-type oscilla- 
tions clearly observable. This is by no means a dangerous approach 
since the maximum deviation from the glide slope, occurring primarily at 
the beginning of the approach, is approximately 17 m. The wavelength of 
the initial wave in the aircraft trajectory measures approximately 2,300 
m (peak to peak), with the remaining waves approximately 1,400 m. The 
phugoid period of the Queen Air (approach speed V a Q = 56 m/s) as deter- 
mined by T = v/27rVa 0 /g from Etkin [9] is 25.5 seconds. The Queen Air 
phugoid wavelength is x = Va 0 T or 1,438 m. However, since this is an 
approximate value and since a pure sine wave at the Queen Air's phugoid 
frequency is not specifically input, the exact Queen Air phugoidal 
wavelength is not expected. The phugoidal-type oscillations are aug- 
mented primarily by the initial increase in head wind from 7.21 m/s to 
8.84 m/s and subsequent decrease to 4.91 m/s. The increasing head wind 
causes the initial rise above the glide slope and an initial increase 
and then decrease of airspeed. Shortly after the airspeed begins to 
decrease, the head wind begins to decrease and the airplane descends 
through the glide slope. The initial encounter with this wave sets off 
the subsequent oscillations where the aircraft descends and picks up 
airspeed, then begins to ascend until it once again loses airspeed and 
begins to descend again. 

Figure 3 shows the trajectory of a simulated Queen Air under auto- 
pilot control descending through the same wind field (Queen Air 17/2). 
Predictably, the trajectory is perfectly aligned with the glide slope 
and all phugoidal-type oscillations are removed. Figures 4 and 5 show 
the minor elevator and thrust control inputs necessary to maintain the 
approach along the glide slope. The sign conventions for elevator angle 
are shown in Figure 6. Since the autopilot represents a flight path 
optimization control system, the initial increasing head wind is 
countered with slightly increased thrust to maintain ground speed and an 
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Figure 4 Autopilot elevator inputs to Queen Air in flight through 
Queen Air flight 17/run 2 wind profile. 
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Figure 5 Autopilot thrust inputs to Queen Air in flight through Queen 
Air flight 17/run 2 wind profile. 



elevator setting moving from negative to positive to pitch the nose down 
to counter the increased rise above the glide slope. As the head wind 
begins its drop from 8.84 m/s to 4.91 m/s, the initially added thrust is 
removed to maintain ground speed and the nose is pitched up to arrest 
descent. The airspeed is still high due to the previously added airspeed 
due to the thrust increase during the increased head wind. 

Throughout the flight, the autopilot tries to closely monitor the 
aircraft's state (airspeed, ground speed, descent rate, and glide slope 
alignment) and retrim the aircraft following the changes in the wind. 

Figures 7, 8, and 9 show trajectory, elevator angle, and thrust 
setting for a simulated Queen Air flown by a pilot with a 70 percent 
performance rating. The control inputs are small and similar to those 
of the autopilot. This is reasonable since the wind shear in this case 
is neither large nor severe. Again the aircraft is flown straight down 
the glide slope with no discernible deviations. Next, a pilot having a 
lower performance rating of 8 percent is examined. Figures 10, 11, and 
12 show the aircraft trajectory and the elevator and thrust control 
settings. The pilot attempts to retrim to the changing wind trend, but 
due to slower response time, the control excursions are larger than 
those of the autopilot and the 70 percent rated pilot. The trajectory 
shows that the pilot cannot adequately counter the initial head wind 
increase, and a small rise above the glide slope can be seen. However, 
the pilot manages to get back to the glide slope, suggesting that the 
main concern in maintaining control is the initially encountered wave 
and that controlling the remainder of the run is relatively straight- 
forward. Finally, a simulated Queen Air commanded by a pilot having a 
5.5 percent performance rating was flown through the wind field. The 
performance rating in this case appears to be the threshold of minimum 
control because, as shown in Figure 13, the pilot is unable to counter the 
effects of the initial wave. The aircraft oscillates about the glide 
slope until finally near the end of the run the pilot is able to bring 
the plane back to intercept the glide slope. The elevator and thrust 
control inputs are presented in Figures 14 and 15, respectively. 
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Figure 7 Simulated trajectory of a Queen Air controlled by a 70 percent 
performance rated pilot flown through Queen Air flight 17/run 2 
wind profile. 
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Figure 8 Elevator inputs by 70 percent rated pilot to Queen Air flown 
through Queen Air flight 17/run 2 wind profile. 
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Figure 11 Elevator inputs by 8 percent rated pilot to Queen Air flown 
through Queen Air flight 17/run 2 wind profile. 
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Thrust inputs by 8 percent rated pilot to Queen Air flown through 
Queen Air flight 17/run 2 wind profile. 
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Figure 13 Simulated trajectory of a Queen Air controlled by a 5.5 percent 
performance rated pilot flown through Queen Air flight 17/run 2 
wind profile. 
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Figure 14 Elevator Inputs by 5.5 percent rated pilot to Queen Air flown 
through Queen Air flight 17/run 2 wind profile. 
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Figure 15 Thrust inputs by 5.5 percent rated pilot to Queen Air flown 
through Queen Air flight 17/run 2 wind profile. 



The results of this simulation show that, although the aircraft is 
disturbed by the wind field in the fixed-stick case, the shear does not 
appear to be strong enough to cause any difficulty to a human-controlled 
vehicle. This is verified by the extremely low performance rated pilots 
who are able to make excellent approaches. 

5.2 Simulated Boeing 727 Fligh t Through Queen Air Flight 17/Run 2 

Wind Field 

An aircraft having the characteristics of a Boeing 727 is flown in 
simulated flight through the same wind field that the Queen Air encoun- 
tered. The time history wind data, which have been converted from range 
data using the Queen Air's average airspeed, were reconverted to range 
data so that both aircraft encounter essentially the same wind field. 

The approach speed used for the B727 runs is 72 m/s as opposed to 56 m/s 
for the Queen Air. 

Figure 16 shows the fixed-stick flight of the simulated B727 through 
the Queen Air wind field. As with the Queen Air, the B727 encounter with 
the initial increasing and decreasing head wind results in a phugoidal- 
type oscillation in the trajectory. These oscillations are of higher 
amplitude and lower damping than those of the Queen Air. This is in 
direct agreement with McCarthy, Blick, and Bensch [3] and Frost and 
Crosby [4], who suggest that large commercial transports will experience 
greater difficulty in flight through wind shear than smaller aircraft 
due to the effects of higher landing speed and higher lift-to-drag ratio 
influence on the exponential phugoid damping period. Table 2 shows 
comparative calculations for the B727 and Queen Air exponential phugoid 
damping periods. The greatest deviation from the glide slope in Figure 
16 is approximately 43 m. The phugoid time period for the B727 given by 
T = /2TrVa 0 /g is 32.5 sec and the wavelength x = V a Q T is 2,341 m as 
compared with the 2,743 m (peak to peak) waves of the simulated trajectory. 

The trajectory of a simulated B727 flown by autopilot through the 
wind field is shown in Figure 17. Predictably the flight is smooth and 
perfectly aligned with the glide slope. In Figure 18 it can be observed 
that the response of the autopilot to the initial head wind increase is 
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TABLE 2. Comparison of Exponential Phugoid Damping Periods (T) of Queen 
Air and Boeing 727. 


T = the time for the phugoid oscillation to decay to 1/e of its initial 
ampl itude. 



where 



Boeing 727 

Queen Air 

V = approach airspeed (m/s) 
0 

72.000 

56.000 

2 

g = gravitational acceleration (m/s ) 

9.800 . 

9.800 

C L = lift coefficient (dimensionless) 
0 

1 .360 

0.639 

C Q = drag coefficient (dimensionless) 
u o 

0.139 

0.080 

T = exponential phugoid damping period (sec) 

71.700 

45.900 
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Simulated trajectory of a B727 under autopilot control flown 
through Queen Air flight 17/run 2 wind profile. 
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Figure 18 Autopilot elevator inputs to B727 in flight through Queen Air 
flight 17/run 2 wind profile. 



a positive deflection in elevator angle. This pitches the nose down to 
counter the nose-up attitude caused by the increased head wind. However, 
in Figure 19 it is noticed that the autopilot optimization control loop 
has decreased thrust when encountering the increased head wind at the 
beginning of the run. With the Queen Air, thrust was increased to 
maintain ground speed. In the case of the B727, the higher approach 
airspeed into the increasing head wind causes too much additional lift 
which tends to force the airplane above the glide slope. Thus, the 
thrust was cut back to prevent the climb above the glide slope. It 
should also be noted that the elevator and thrust control inputs to the 
B727 are of larger magnitude than those of the Queen Air. 

The trajectory of an aircraft controlled by a pilot with a 70 
percent performance rating is given in Figure 20. The aircraft flies 
straight down the glide slope. Figures 21 and 22 provide the elevator 
and thrust control inputs which are similar in magnitude to those made 
by the autopilot. 

The flight of an aircraft flown by a pilot with a low performance 
rating (5.5 percent) is next considered. The flight profile (Figure 23) 
shows that this pilot is unable to fully counter the initial increased 
head wind and the aircraft rises slightly above the glide slope at the 
beginning of the run. Following a slight descent below the glide slope 
due to the decreasing head wind, the pilot manages to correct and main- 
tain a near-perfect approach. Note that the elevator and thrust controls 
(Figures 24 and 25) show less "fine tuning" than those of the 70 percent 
rated pilot. 

This above-described simulation shows that although the aircraft 
oscillates due to the wind shear for the fixed-stick run, there is no 
problem in maintaining an excellent approach in a piloted aircraft even 
for a pilot with a low performance rating (5.5 percent). 

5.3 Simulated B727 Flight Through Sinusoidal 6 m/s Amplitude Wave Near 

Phugoid Frequency 

A simulation is carried out flying the simulated B727 through a 6 
m/s head wind shearing to a 6 m/s tail wind. This wave disturbance is 
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Figure 21 Elevator inputs by 70 percent rated pilot to B727 flown through 
Queen Air flight 17/run 2 wind profile. 
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Figure 22 Thrust inputs by 70 percent rated pilot to B727 flown 
Queen Air flight 17/run 2 wind profile. 
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Figure 24 Elevator inputs by 5.5 percent rated pilot to B727 flown through 
Queen Air flight 17/run 2 wind profile. 
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Figure 25 Thrust inputs by 5.5 percent rated pilot to B727 flown through 
Queen Air flight 17/run 2 wind profile. 



similar to that of the Queen Air 17/2 case except that a wave at approxi- 
mately the B727's phugoid frequency is encountered by the aircraft. 

Such a wave disturbance is used to simulate the thunderstorm cell wind 
environment which has been the cause of either an accident or, at 
minimum, an incident to several aircraft which encountered similar 
waves. An aircraft flying through a thunderstorm initially flies into a 
head wind (the thunderstorm cell outflow). This longitudinal head wind 
component diminishes to zero as the aircraft flies through the central 
down draft region of the cell and becomes a tail wind as the aircraft 
flies into the outflow on the other side. 

Figure 26 shows the fixed-stick trajectory of the simulated B727 
flying into a constant 6 m/s head wind which begins to decrease at 
X/HA = 8.3 to become a 6 m/s tail wind at X/HA = 21.2. The aircraft 
drops below the glide slope and loses airspeed until X/HA = 14.8 when 
the airspeed starts to increase again. However, since the head wind is 
still decreasing, the aircraft continues to fall as the airspeed continues 
to build up. Finally, the aircraft then begins to climb at X/HA = 33. 

The encountered wave disturbance thus sets off phugoidal oscillations in 
the trajectory. Since the sinusoidal wind velocities are computed using 
a time scale based on aircraft position along the glide slope and average 
approach speed, the sine wave is approximate since the aircraft does not 
travel collocationally along the glide slope during the wave encounter. 

The flight path of a simulated B727 flown by the autopilot is given 
in Figure 27. The beginning of a very slight descent below the glide 
slope occurs at X/HA = 11, which is the same location that the fixed-stick 
aircraft (Figure 26) begins to drop. The autopilot, however, quickly 
begins to pitch the aircraft nose-up (Figure 28) and increases the thrust 
(Figure 29) to counter the loss of airspeed and deviation below the 
desired glide path. The aircraft comes back up to intercept the glide 
slope and continues the approach through the remaining 6 m/s tail wind 
without deviation. Note that the retrim for maintaining flight through 
the tail wind requires a higher thrust setting than for the initial head 
wind as well as an elevator angle setting of approximately -8°. Once 
again it can be observed that by correct countering of the wave-form wind 
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disturbance, the oscillations of the fixed-stick case can be almost com- 
pletely eliminated. 

Simulated landings with pilots of various performance ratings were 
flown down the same sinusoidal wave. Figure 30 illustrates that a pilot 
with a 25 percent rating can easily negotiate the disturbance with only 
a slight dip below the glide slope. The wave-like "humps" in the elevator 
control setting (Figure 31), observed from X/HA = 0 to 4 and from X/HA = 

35 to 80, represent pilot inputs to counter the aircraft's slight natural 
phugoidal oscillation in the steady 6 m/s head wind before the wave and in 
the tail wind after the wave, respectively. 

Figure 32 shows the flight path of an aircraft controlled by a 
simulated pilot with a 5.5 percent performance rating. Passing through 
the shear wave, this pilot does not react quickly enough to the decreasing 
airspeed caused by the decreasing head wind, and the aircraft falls a 
maximum of 23 m below the glide slope. The pilot then regains control and 
applies step inputs to the elevator and throttle (Figures 33 and 34) to 
get back on the glide slope. In general, the pilot's approach is fairly 
good. 

The results of this simulation show that a 6 m/s thunderstorm 
environment phugoid wave is not very severe from the standpoint of pilot 
control, as even a pilot with low performance capability can make the 
approach without too much deviation from the glide slope. However, had 
the deviation below the glide slope occurred relatively close to the 
runway threshold, the 5.5 percent rated pilot would have been forced to 
execute a go-around. 

5.4 S imulated B727 Flight Through Sinusoidal 14 m/s Amplitude Wave Near 

Phugoid Frequency 

In the following simulations, the B727 is flown through a longitudi- 
nal wave disturbance representative of a very strong thunderstorm cell 
wind profile. A 14 m/s amplitude wave is chosen based on longitudinal 
thunderstorm wind data analyzed by Frost, Camp, and Wang [10]. The air- 
craft is initially trimmed for flight in a steady 14 m/s head wind which 
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Figure 30 Simulated trajectory of a B727 controlled by a 25 percent rated 
pilot flown through a 6 m/s amplitude phugoidal sine wave. 
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Figure 31 Elevator inputs by 25 percent rated pilot to B727 flown through a 
6 m/s amplitude phugoidal sine wave. 
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Figure 33 Elevator inputs by 5.5 percent rated pilot to B727 flown through a 
6 in/s amplitude phugoidal sine wave. 






shears as a wave at approximately the aircraft's phugoid frequency to a 
steady 14 m/s tail wind for the remainder of the run. 

The simulated trajectory of an autopilot-controlled aircraft is shown 
in Figure 35. The aircraft encounters the decreasing head wind at approxi- 
mately X/HA = 8 and begins to fall below the glide slope. The autopilot 
increases thrust (Figure 36) and moves the elevator to full negative 
deflection (Figure 37) to pitch the aircraft nose-up. The aircraft falls 
a maximum of 15 m below the glide slope until it responds to the added 
controls and pulls out of its descent at X/HA = 14. Due to the added 
controls, the aircraft ascends above the glide slope to a maximum of 23 m 
at X/HA = 30. The autopilot then decreases thrust to allow the aircraft 
to gradually fall back and recapture the glide slope and finish with a 
near perfect approach. While it is unlikely that an aircraft would 
execute an approach in a 14 m/s tail wind, note that the autopilot trims 
the aircraft with full nose-up elevator deflection and twice the initial 
thrust to maintain flight in this tail wind. 

Figure 38 shows the simulated trajectory of an aircraft controlled by 
a pilot with a performance rating of 50 percent. The aircraft drops 
approximately 18 m below the glide slope when it encounters the decreasing 
head wind. Again, thrust is increased (Figure 39) and elevator angle is 
increased to full nose-up deflection (Figure 40) to counter the effect of 
the decreasing airspeed. However, due to the farther drop below the glide 
slope and slower response time of the pilot as compared with the auto- 
pilot, the aircraft does not rise as high above the glide slope as it did 
under autopilot control. Thrust is then decreased so that the aircraft 
can descend to intercept the glide slope and maintain an excellent 
approach through the remaining tail wind. 

The flight profile and control inputs of an aircraft controlled by a 
25 percent rated pilot (Figures 41 through 43) are examined next. During 
the encounter with the shear wave, the aircraft falls a maximum of 23 m 
below the glide slope. This would obviously create a hazardous situation 
had the aircraft encountered such a wave in the immediate vicinity of the 
runway threshold. Because of the large drop below the glide slope and 
slow response of the pilot, the aircraft rises only slightly above the 
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Figure 35 Simulated trajectory of a B727 under autopilot control flown 
through a 14 m/s amplitude phugoidal sine wave. 
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Figure 36 Autopilot thrust inputs to B727 in flight through a 14 m/s 
amplitude phugoidal sine wave. 
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Figure 37 Autopilot elevator inputs to B727 in flight through a 14 m/s 
amplitude phugoidal sine wave. 
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Figure 40 Elevator inputs by a 50 percent rated pilot to B727 flown 
through a 14 m/s amplitude phugoidal sine wave. 
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Figure 41 Simulated trajectory of a B727 controlled by a 25 percent rated 
pilot flown through a 14 m/s amplitude phugoidal sine wave. 
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glide slope as it pulls out of its descent. Thrust is decreased at X/HA 
= 21 to return the aircraft to the glide slope where it remains for the 
length of the run. 

Finally, the performance of a 5 percent rated pilot in an aircraft 
encountering the 14 m/s shear wave is examined. Figures 44, 45, and 46 
show the flight profile and the elevator and thrust inputs, respectively. 
As the aircraft enters the decreasing head wind, the airspeed decreases 
and the aircraft begins to descend. The pilot's response is extremely 
slow, allowing a very high sink rate to build up. Because of the large 
lag on the controls and severe change in the wind, control of the aircraft 
is lost. The initial loss of control builds up to the point where the 
control loop starts to decrease thrust between X/HA = 13 to X/HA =17 
(Figure 46). At this point there appears to be a small "bug" in the control 
loop. The thrust should not decrease with increasing sink rate even 
though the aircraft is experiencing an increasing speed. 

It is believed that a small correction to the gain on the vertical 
velocity deviation in the pitch control variable e c , which has been kept 
constant for all previous cases, will correct this slight anomaly. 

However, the severity of this wind field and the low performance rating 
used in this case is apparently responsible for this effect. The proper 
setting of this gain will require further investigation. Nevertheless, 
the aircraft does eventually begin to regain thrust control and pull out 
of its descent at approximately X/HA = 18 when the vertical velocity 
peaks, compensating for the low gain. Thrust is eventually increased to 
a higher setting (Figure 46) and step inputs are made to both elevator and 
thrust to return the aircraft to the glide slope. The response of the 
aircraft to the step inputs is shown by the wave-like trajectory in 
Figure 44 from X/HA = 20 to X/HA = 75. 

The results of this simulation illustrate the severity of strong 
thunderstorm-related shear in the horizontal wind. The autopiloted and 50 
percent rated pilot-controlled aircraft maintained fairly good approaches 
but not without significant deviation below the glide slope. The 25 
percent rated pilot-controlled aircraft experienced a fairly large drop 
below the glide slope of such magnitude that a go-around would have had to 
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Figure 45 Elevator inputs by a 5 percent rated pilot to B727 flown 
through a 14 m/s amplitude phugoidal sine wave. 
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Figure 46 Thrust inputs by a 5 percent rated pilot to B727 flown through 
a 14 m/s amplitude phugoidal sine wave. 




be executed had the shear wave been encountered relatively close to the 
runway threshold. Finally, the inability of the low performance 5 percent 
rated pilot to negotiate the high rate of descent caused by the shear wave 
shows the danger of slow response and failure to properly recognize the 
effects of wind shear on approach. 

5.5 Simulated B727 Flight Through Full Sinusoidal 10 m/s Amplitude Wave 

Near Phugoid Frequency 

In the final simulations the B727 is flown through a full 10 m/s 
head wind to tail wind sinusoidal phugoid frequency wave. The wind 
increases from 0 m/s to become a 10 m/s head wind which shears to a 
10 m/s tail wind and then returns to zero wind for the remainder of the 
run. 

Figure 47 shows the trajectory of a simulated B727 under autopilot 
control flown through the sinusoidal wind field. The aircraft is forced 
slightly above the glide slope when it encounters the increasing head 
wind and drops slightly below the glide slope as it enters the increasing 
tail wind. The aircraft then rises above the glide slope as the tail 
wind diminishes to 0 m/s. The autopilot eventually returns the aircraft 
to the glide slope to land on target. 

Figures 48, 49, and 50 show the simulated flight path, elevator 
inputs, and thrust control inputs, respectively, for a piloted B727 
flown through the phugoidal sine wave. The pilot model used for this 
run was based on the data for pilot A from Adams and Bergeron [7] 

(Table 1). This pilot was rated relatively high based on high static 
gains, high lead time constants, and low lag time constants. The 
aircraft rises above the glide slope when encountering the increasing 
head wind. The pilot cuts thrust to zero and pitches the nose down to 
counteract the rise above the glide slope. The aircraft then encounters 
the head wind to tail wind shear and drops below the glide slope, due 
to the previously decreased power, necessitating a thrust increase and 
pitch-up of the nose. Eventually, the aircraft is brought back to the 
glide slope and retrimmed for the zero wind condition for a safe landing. 
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Figure 49 Elevator inputs by pilot A to B727 flown through 10 m/s amplitude 
phugoidal sine wave. 
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Figure 50 Thrust inputs by pilot A to B727 flown through 10 m/s amplitude 
phugoidal sine wave. 



Figure 51 shows the trajectory of a simulated B727 flown through the 
shear wave by pilot F. As shown from the data of Adams and Bergeron [7] 
in Table 1, this pilot is of relatively low skill due to a low static 
gain, low lead time constant, and high lag time constant. The aircraft 
as controlled by this pilot rises highest above the glide slope and 
drops farthest below the glide slope due to higher response lag time. 

The pilot is unable to bring the aircraft back to the glide slope and 
lands short. The wave- like humps in the trajectory represent attempts 
to stabilize the approach angle. 

This simulation compares the relative degrees of control applied to 
an aircraft landing through a phugoidal shear wave. The optimization 
autopilot rapidly corrects for the wind changes with low lag time 
and maintains the tightest control. Pilot A possesses higher lag than 
the autopilot and illustrates the danger of decreasing thrust when 
encountering a head wind to tail wind shear. The aircraft dropped below 
the glide slope but the pilot was able to regain control and recover the 
approach. The lower skilled, higher' lagged pilot F was unable to safely 
counter the shear and landed short. 
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6.0 CONCLUSIONS 


The purpose of this study was to evaluate the response of the 
pilot-aircraft system to wind shear. This was carried out by performing 
a parametric study of pilot performance rating and by use of actual 
pilot data. Performance ratings ranged between 100 percent autopilot 
control (perfect pilot) and "hands off" or 0 percent control (fixed-stick), 
with the human pilot having any performance rating in this range. Fixed- 
stick, autopilot, and manned performance of both a simulated B727, a 
medium-sized commercial transport, and a simulated Queen Air, a small 
commuter-type aircraft, flown through a glide slope longitudinal wind 
profile detected by Doppler radar were compared. The wave-form wind 
disturbance was shown to excite the phugoid oscillations of both air- 
craft when flown in the fixed-stick mode but presented no control problems 
for manned aircraft. The fixed-stick simulations also showed agreement 
with the theories of McCarthy, Blick, and Bensch [3] and Frost and Crosby 
[4] stating that the oscillatory effects of a shear wave are greater on 
transport aircraft than on smaller aircraft as predicted by aircraft 
stability theory. 

A simulation was also carried out flying the fixed-stick, autopilot 
controlled, and manned aircraft with characteristics of a B727 through a 
6 m/s amplitude head wind to tail wind phugoid frequency shear wave 
representative of the thunderstorm cell wind environment. This case also 
revealed phugoid oscillations but clearly illustrated that this shear 
wave was not a serious problem to a manned vehicle. 

However, in the stronger disturbance of a 14 m/s head wind to tail 
wind phugoid frequency shear wave, significant deviation from the glide 
slope was noted for the autopilot, 50 percent rated pilot, and 25 percent 
rated pilot flight simulations. The low performance, 5 percent rated pilot 
initially lost control of the aircraft and dropped farthest below the 
glide slope. However, thrust was eventually increased to bring the 
aircraft back up to the glide slope. 
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In flight simulations through a full 10 m/s phugoidal sine wave, 
comparisons were made between autopilot control and control by simulated 
pilots of varying skill. The autopiloted aircraft executed the best 
approach, while the high-skilled pilot descended below the glide slope but 
was eventually able to bring the aircraft back onto the glide path. 
However, the low-skilled pilot could not maintain adequate control and 
landed short. 

The four wind profiles utilized in this study do not by any means 
limit the severity or variety of conditions to be found in the atmosphere. 
However, this research provides a better understanding of pilot-aircraft 
response to longitudinal wind shear and serves as an important input 
toward the development of a Doppler radar wind shear warning and 
detection system. 
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APPENDIX A 


DERIVATION OF THRUST AND ELEVATOR ANGLE 
CONTROL FUNCTIONS ( Fy , 6 E ) 

The three degrees of freedom of the two-dimensional aircraft motion 
model are given by the following equations from Frost and Crosby [2] 
q = D 7 F t ♦ D s V^C m (1) 

yV = Dg sin ( 6y + a)Fy - D 2 cos y + D.jV 2 (C l cos 6 - Cp sin 6) (2) 

V = Dg cos (<5y + a)Fj - D 2 sin y - D-j V 2 (C E sin 6 + Cp cos <S) (3) 

Equations 1, 2, and 3 represent pitch motion about the y axis, vertical 
motion along the z axis, and longitudinal motion along the x axis, 
respectively. The axis system has the origin fixed at the airplane's 
center of mass with the x axis oriented along the ground speed vector. 
Expanding C m , C^, and C D plus adding an arbitrary pitch control variable, 
0g , to Equation 2 yields 

q = D 7 F t + DrV?(C + C a' + C + (D,C /Vjq + (D.C /Vja‘) (4) 
^ 7 T 5 a rn m m~ E '4m' a ^ ' 4 m> a ' 

O a q a 

y V = Dg sin (<5y + a)Fj - D 2 COS y + D i V^ C C|^ + C^ a 1 + C^ <5^ 

0 a 

+ (D 4 C l /V a )q + (D 4 C L ./V a )a') cos 5 - D^ 2 ^ + Cp a 1 
q a O a 

+ C n a 12 ) sin 6 - e r (5) 

U o L 

a z 

V = Dg COS (6y + a)Fy - D 2 sin y - D-|V 2 (C E + C^ a 1 + Cy 6^ 

0 a 6 ^ 

+ (D 4 C l /V a )q + (D 4 C L ./V a )a') Sin 6 - D^ 2 ^ + Cp a 1 + C Q a' 2 ) 
q a 0 a a 2 


x COS 6 
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( 6 ) 



Collecting coefficents of Fy, a 1 , and in the three equations and 
setting all time derivatives to zero 

K V a D 5 c m 

a 

Dg sin (5 t + a) - V2 (D-,C d + 2 a') sin 6 + D-,C L V 2 cos 6 

D 6 cos (6 t + a) - V 2 ( *') cos 5 - V 2 sin 6 


5am, 


D 1 V 2 C, cos <5 a 1 

6 e 

- D 1 V I C I si " 6 6 e 
6 E J l J 

Next, separating the remaining terms of the three equations, gives 


■-< D 5 C m V a + WrO + D 5 D 4 C m.“'> 
0 q a 


Kff'c} 


~( D 1 C L V a + D 1 D 4 C L q + D 1 D 4 C L.“^ cos 6 + D 1 C D V a sin 6 
L 0 q a 0 



I 


For computational convenience and clarity, we write Equations 7 and 8 as 



To solve for Fy and , we take the inverse of the 3x3 matrix on the left- 
hand side of Equation 9 which is 

'< G 5 G 3 - W V a - (G 8 G 3 - G 2 G 9 )V a < G 8 G 6 ' G 5 G 9> V 5 ' 

SET -< G 7 G 3 - G 4 G 6^a (G 1 G 3 ’ G 4 G 9> V a - (G 1 G 6 ’ G 7 G 9 )V a < ,0 > 

< G 7 G 2 - G 4 G 5> V a -< G 1 G 2 ‘ G 4 G 8 )V a (G 1 G 5 ‘ G 7 G 8> V a 

where the determinant, DET, is equal to 

(G1G5G3 + G g GgG 4 + G g G 7 G 2 - G 4 G 5 G g - G 2 G 6 G 1 “ 6 3 G 7 G 8 ^ V a 

[ Q 2 (Gg G 7 " G 6 G 1^ + G 3^ G 1 G 5 “ G 7 G 8^ + G 4^ G 6 G 8 " G 5 G 9^ V a 
= -(s 2 h, - G 3 H 2 + G 4 H 3 )V^ 

■ -(«4> V a 
* - (H 5>^ 


Multiplying both sides of Equation 9 by Equation 10 yields 



Since we are solving for Fy and 5^, the terms of the second row of the 
inverse corresponding to a 1 were not derived. Solving for Fy, 
and <5g gives 
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(12) 


F I ' K Tl V a + K T2 V + K T3 V + 


6 E = K Dl V a + K D2 V + K D3 V + 


K T4 6 C 
K D4 9 C 


where 


K T i 

1 

" H n G i2 + H io G n + H 9 G io 

K T2 


-h 9 d 2 

K T3 


- H 10 D 2 

K T4, 


“ H 10 


(13) 


(14) 


and 


K D1 

l 

H 8 G 1 2 " H 7 G 1 1 " H 6 G 10 

K D2 


Hg D 2 

K D3 


H 7 D 2 

K D4 


H 7 


(15) 


Equations for G-, H., D. , and C appear in detail in the following 
T j k n 

section. 


Equation Form of G.., H ^ , D k> and C n 
G 1 = °7 

G 2 = -(C-j + C 2 a') cos 6 - C 3 sin 6 
G 3 = -C 4 sin <5 
G^ = Dg cos (5j + a) 

Gg = Cg COS 6 - (C-j + C 2 a') Sin 5 

Gg = cos <5 

G 7 = Dg sin (Sy + a) 
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G 8 C 5 
G 9 = C 6 

G 10 = C 7 V a cos 6 + ( C 8 V a + C 9 q + C 10“^ sin 6 
G 11 = (C 8 V a + C 9 q + C 10“' ) cos 6 “ C 7 V a sin 6 
G 1 2 = C ll V a + C 12 q + C 13 a ' 

H ] = G 6 G 1 " G 7 G 9 

H 2 = G 5 G 1 _ G 7 G 8 

H 3 = G 5 G 9 “ G 6 G 8 

H 4 = G 2 H 1 " G 3 H 2 + G 4 H 3 

H 5 =V a H 4 

h 6 =h 2 /h 5 

H7 = (G 2 G-| “ 

H g = (G 2 G 7 - G 4 G 5 )/H 5 

H9 = 

H 10 = (G 2 G 9 " G 3 G 8 )/H 4 

H 11 = ^ G 2 G 6 “ G 3 G 5^ H 4 


D 


= (p/2) (SHA/m) 

D 2 = (gHA/V2) 

D 

D 4 = (C/2HA) 

D 

D 6 = (gHA/mV 2 ) 

°7 = < l tW 

The following are the "C" coefficients used in variable gain computations 
of Fy and 6^: 


(HA 2 /V2) 


= (p/2) (ScHA 2 /g I yy ) 
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I III 



D 1 C D 2 


D 1 D 4 C L. 


Wl- 

a 


D.D.C 
5 4 m. 


Dj-D-C 
5 4 m* i 

a 
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APPENDIX B 


CHARACTERISTICS OF AIRCRAFT USED IN 
NUMERICAL FLIGHT SIMULATIONS* 



Boeing 727 

Queen Air 

Reference (scale) altitude HA (m) 

91.4 

91.4 

Initial approach airspeed V= (m/s) 

71.9 

56.4 

Initial flight path angle y (deg) 

-3.0 

-3.0 

Aircraft mass (kg) 

63,945.6 

3,469.2 

2 

Moment of inertia I (kg-m ) 

6.1 x 10 6 

7.8 x 10 3 

Moment arm of thrust vector Ly(m) 

0.0 

0.0 

Angle between Fy and fuselage reference 
line (deg) 

0.0 

0.0 

Chord length c” (m) 

5.0 

1.8 

2 

Wing area S (m ) 

145.0 

27.3 


Boeing 727 

Queen Air** 

Cl 0 

1.360 

0.639 

C L 

a 

5.04/rad 

5. 28/ rad 

\ 

0.007/deg 

0.007/deg 


9.3/rad 

2. 9/ rad 


*Courtesy of E. F. Blick, Professor of Aerospace Engineering, 
Department of Engineering, University of Oklahoma, Norman, Oklahoma. 

**Based on approach speed of 62 m/s. 

86 



C L. 

a 

6.6 

1.08 

Cd o 

0.139 

0.08 

C D 

a 

1 .245/rad 

0.33/rad 

C D 2 

0.0/rad 2 

0.0/rad^ 

C 

m o 

0.0 

0.0 

C 

m 

a 

-1 .47/rad 

-1 .0/rad 

c m 

6 E 

-0.025/deg 

-0.025/deg 

c 

m q 

-29.5/rad 

-8.7/rad 

C m- 

a 

-1 .77/rad 

-3.24/rad 
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